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Abstract 



Previous studies of gravitationally bound states of ultracold neu- 
trons showed the quantization of energy levels, and confirmed quan- 
tum mechanical predictions for the average size of the two lowest en- 
ergy states wave functions. Improvements in position-like measure- 
ments can increase the accuracy by an order of magnitude only. We 
therefore develop another approach, consisting in accurate measure- 
ments of the energy levels. The GRANIT experiment is devoted to 
the study of resonant transitions between quantum states induced by 
an oscillating perturbation. 

According to Heisenberg's uncertainty relations, the accuracy of 
measurement of the energy levels is limited by the time available to 
perform the transitions. Thus, trapping quantum states will be neces- 
sary, and each source of losses has to be controlled in order to maximize 
the lifetime of the states. We discuss the general principles of tran- 
sitions between quantum states, and consider the main systematical 
losses of neutrons in a trap. 

1 Introduction 

As predicted by quantum mechanics, a neutron bouncing above a mirror 
has discrete energy states. Due to the extreme weakness of gravity, these 
quantum states exhibit outstanding properties. 

First, the mean height of the n th state is much larger than atomic size, 

z n « (n - l/4) 2/3 x 11 /im, (1) 

which has been used to prove the quantization of energy at the ILL high 
flux reactor in Grenoble [HE]- A neutron flux has been measured through a 
narrow slit between a horizontal mirror and a scatterer above. We observed 
the discrete behaviour of this flux as a function of the height of the slit. In 
particular, the flux is zero if the height is less than 10 /im, since no quantum 
state can penetrate throught the slit. The classical turning points of the two 
first states have then been determined to be [3], 

zT P = 12.2 ± 1.8 sys ± 0.7 stat /im, (2) 

zf = \{l\z\l) = 13.7 /im, 

z? p = 21.6±2.2 sys ±0.7 stat /im, (3) 
3 

4 h = - (2\z\2) = 24.0 /im. 
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Higher states are more difficult to resolve. We expect to be able to improve 
the precision of position-like observables by one order of magnitude at most. 

One second important property of the quantum states is the extreme 
smallness of the energies. Indeed, the energy of the n th state is, within the 
Bohr-Sommerfeld approximation, 

E n « {n - l/4) 2/3 x 1.7 peV. (4) 

Then, the corresponding frequencies f n — ^ are small as well. They are 
in the kilohertz range, which makes it possible to probe the quantum states 
with an oscillating perturbation at these easily accessible frequencies. In the 
next section, we will discuss the principle of the GRANIT experiment, which 
aims to induce resonant transitions between quantum states. The accuracy 
of the transition frequencies is limited by the lifetime in the quantum levels. 
To increase the duration of the perturbation, a trap has to be built, and the 
third section is devoted to the estimation of lifetimes of quantum states in 
this trap. 



2 Resonant transitions in the GRANIT ex- 
periment 

Let us assume that a given neutron bouncing above a horizontal mirror stands 
in a pure initial state \N) concerning its vertical motion. We apply a periodic 
perturbation, 

V(t) = Re (v(z)e iuJt ) , (5) 



induced by an oscillating magnetic gradient, oscillations of the horizontal 
mirror itself, or by the motion of neighbouring mass. If the angular frequency 
uj is close to £ "~ £:jv , which corresponds to the transition to an excited state 
\n), then a Rabi resonance is expected. The probability to observe this 
transition is well known, 

PN^n(t) = — ^ -j sin 2 U (u - u Nn ) 2 + Q 2 Nn I) , (6) 



Nn 



where Qnu = j;{ n \V( z ) \N) is the Rabi frequency, defining the perturba- 
tion strength for the N — > n transition. The maximum of the transition 
probability, 

pmax _ \ 

r N~*n — , s 2 \' ) 
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TRANSITION SPECTRUM FROM STATES 1 AND 2 




f[Hz] 



Figure 1: Maximum transition probability for neutrons prepared in the first 
quantum state (straight line) and for neutrons prepared in the second quan- 
tum state (dotted line), as a function of the perturbation frequency. The 
Rabi angular frequency Q satisfies I = 0.1 s. 



has a Lorentzian shape, and is reached if the perturbation is applied for the 
so-called pulse time, 



J- pulse , V u 7 

yj(u ~ UJ Nn ) 2 + tt 2 Nn 

Fig. [U shows the maximum transition probability of the first two quantum 
states for frequencies below 1 kHz. 

Obviously, the pulse time must be smaller than the storage time of neu- 
trons in a given quantum state. Thus, we can deduce from the previous 
discussion the two fundamental reasons to increase the storage time. First, 
the resonance formula (EJ) contains the Heisenberg's uncertainty relation, 
AE ■ T > h/2, where T is the storage time (T > T pu i se ) and AE the width of 
the resonance curve. As a work hypothesis we consider that the accuracy of 
GRANIT to measure energy transitions equals this resonance width, which 
is then inversely proportional to the storage time. In expression ((HI) we can 
see the second reason to increase the storage time, since the perturbation 
strength needed for a 100% probability at the resonance is inversely propor- 
tional to the pulse time as well. In the remaining of this section we will 
explicit this relation connecting strength of perturbation and storage time in 
the case of magnetically induced transitions, and in the next section, we will 
estimate the lifetime of the quantum states in the GRANIT trap, which is 
of primary importance. 

For now, let us list the characteristic timescales of the problem: 
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1. The pulse time needed to resolve the states is about 10 ms, leading to a 
resonance curve with a width of half the separation frequency between 
neighbouring states. 

2. The flow through mode time is about 75 ms, which is the time a neutron 
of horizontal velocity 4 m/s takes to pass above a 30 cm long horizontal 
mirror. Of course this time could be increased with slower neutrons or a 
longer mirror, but we need to trap neutrons to gain orders of magnitude. 
Let us notice, however, that it is possible to resolve the states in flow 
through mode. 

3. The ultimate storage time is given by the j3 decay lifetime of the neu- 
tron, i.e. 886 s. If neutrons can be trapped in a given quantum level 
that long, energy levels could be measured with a relative accuracy of 
about 10 -6 . 

4. We may consider the possibility of radiative decays of the quantum 
states by spontaneous graviton emissions. Nonetheless the correspond- 
ing time has been found to exceed the age of the universe by several 
orders of magnitude [I], and can be completely neglected. Therefore, 
contrary to quantum levels in atomic or nuclear physics, the gravita- 
tional quantum levels are absolutely stable, except for external pertur- 
bation or the decay of the neutron itself. 



2.1 Transitions induced by magnetic field 

Let us examine how strong this perturbation has to be in the case of magnet- 
ically induced transitions. Since an uniform magnetic field does not couple 
different quantum states, we apply the magnetic gradient 

B = (/3 z e z + f3 x e x )z cos(cjt) , (9) 

whose variation with time induces two different transitions. The first one 
does not change the spin-state and is induced by the gradient of the vertical 
component, 

V(t) = -fizPz z cos(wt), (10) 

where /i is the magnetic moment of the neutron. Contrary, the second tran- 
sition is induced by the gradient of the horizontal component and changes 
the spin-state of the neutron, 

Vaip(t) = -pvPx z cos(ut). (11) 

1 This field actually violates Maxwell's equations. We may add the gradient term 
{flx^z — fte x ) xcos{u>t), but this term does not couple different quantum states. 
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Figure 2: Magnetic field gradient needed for a 100% transition probability at 
resonance as a function of the time during which the transition is performed, 
for the 1 —>■ 2 and the 2^7 transition. Also shown is the minimum time 
needed to resolve the quantum states, the ultimate j3 decay time, and the 
typical time in flow through mode. 



This property of spin-flipping is important since it can be used to detect that 
a transition has occured. In any case, the Rabi angular frequency is given by 

^7T' lC = \^ (n\z\N). (12) 

Fig. [5] shows the magnetic gradient needed for a 100% transition probability 
at the resonance as a function of the pulse time. If a storage time of 10 
s is achieved, then we can use a magnetic gradient as low as 0.01 T/m to 
induce the 2^7 transition. In flow through mode, a 100% probability for 
the 1 — >• 2 transition can be induced by a magnetic gradient of 0.1 T/m. By 
using the spin-flip, the sensitivity to detect the transition can be increased 
significantly permitting to observe transitions even with smaller gradients. 

2.2 Effects of the Earth' rotation 

Besides the (3 decay lifetime of the neutron, there is another (almost) un- 
avoidable effect limiting the precision of the transitions energies. Indeed, the 
rotation of the Earth induces non-inertial effects, described by the potential 
— Osarth • L, where OEarth is the rotation vector of the Earth, and L is the 
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angular momentum of the neutron relatively to the Earth' center. The main 
effect induced by this potential is the coupling between horizontal velocity 
and vertical motion, 

^Earth rotation — ^EarthCOS# TUV^s % (13) 

where vns is the neutron velocity in North-South direction and 9 is the 
latitude (in Grenoble cos# = 0.7). Thus, the effective vertical acceleration 
of the neutrons is shifted, depending on its horizontal velocity. For the 
first energy state, we can estimate the shift between neutrons travelling in 
North-South direction and neutrons travelling in East- West direction, for an 
absolute velocity of 5 m/s, 

m) „ io-«. (u) 

\ &1 / Earth rotation 

Figure E] shows the corresponding frequency shifts as a function of the state 
number. This effect is close to the ultimate sensitivity for energy levels 
measurements for the ground state, and is even higher for excited states. 
Let us notice that Earth rotation also induces a Zeeman shift between spin- 
up and spin-down states, but this effect is much smaller than the ultimate 
sensitivity, 

(A£) Earth zeeman = ^Earth « 6 X 10~ 8 P eV. (15) 



3 Lifetimes of quantum states in the trap 

The trap of quantum states in the GRANIT setup will look like in fig. [31 
It is a 30 cm square bottom horizontal mirror surrounded by vertical side 
walls. The horizontal velocity of neutrons in the trap will be about 5 m/s, 
this velocity is limited by the Fermi potential of the side walls. In this sec- 
tion, we estimate the loss rate of neutrons due to geometrical imperfections 
of the trap as shown in fig. [3j that is, the waviness of the bottom mirror, 
the deviation from verticality of the side walls, and the corner defects. Other 
sources of losses, such as the seismic noise, a possible remaining inhomoge- 
neous static magnetic field, interaction with magnetic impurities under the 
surface, interaction with dust or a hydrogen layer on the surface, interaction 
with low-energy phonons in the mirror and diffuse scattering, will not be 
considered here. 



3.1 Losses due to waviness of the bottom mirror 

Let us assume that the bottom mirror has a wavy profile £(x), so that a 
neutron with horizontal velocity v sees in its rest frame a time-dependant 
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DEVIATION 
FROM 

VERTICALITY 



L = 0.3 m 

Figure 3: Scheme of the GRANIT trap for gravitational quantum levels 
of UCN. It shows the main geometrical sources for losses of neutrons: the 
waviness of the bottom mirror, the deviation from verticality of the side wall, 
and the corner defect. 



boundary £(vt). We aim to calculate the loss rate due to the transfer of 
horizontal velocity to vertical motion. This problem has already been solved 
in the case of seismic noise [5], and will be developed in more details in a 
forthcoming publication. Following this perturbation approach, it can be 
shown that the rate for the iV — > n loss channel reads 



pwavy _ ( m 9\ 2 1 pQn ( f. 



Nn 



(16) 



where f^n is the frequency of the N — > n transition and PSD(K) is the 
power spectral density of spectral noise, 



PSD(K) = lim 



( £(x)e 2i7TKx dx 
Jo 



(17) 



The GRANIT bottom mirror has not been nuilt yet, but the power spectral 
density of a high-quality 300 mm Si substrate has been measured with several 
characterisation methods [6], 

/ K \ ~ 2 - 9 

PSD(K) = 2 x 10~ 4 nm 2 mm. (18) 

VI mm -1 / 
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It is thus possible to estimate the loss rate of a given quantum level, summing 
all final states in eq. ffTUj) . The result is shown in fig. [Ul as a function of 
initial state quantum number, and for a horizontal velocity of 5 m/s. For the 
first 30 quantum states of interest, the rate of losses due to the waviness of 
the bottom mirror is much smaller than the j3 decay rate. 

3.2 Losses due to deviation from verticality of the side 
wall 

A vertical side wall is obviously a necessary feature of the GRANIT trap. 
Here we estimate how precisely vertical this wall has to be. The probability 
of transition to different states due to deviation from verticality of the side 
walls can be calculated in the framework of the " sudden kick" approximation. 
Indeed, the collision with a wall can be treated as a sudden kick, during 
which the neutron gets the vertical momentum k sin(2a), where a is an angle 
between the wall and the vertical direction, k = mv/h is the wave number 
of horizontal motion. Following the solution given by A.B. Migdal (see e.g. 
[7]), we immediately get the probability to leave a given state \N) during a 
collision, 

P» a11 = 1 - | (N\ exp(ik sm(2a)z) \N) | 2 (19) 

The result is shown in fig. HI for neutrons with horizontal velocity of 5 m/s, 
and for different values of a. Notice that this probability is a quadratic 
function of v a (for small a). As shown in fig. [BJ the corresponding loss 
rate r™ 11 = ^-P™ 11 (for a = 10 -5 rad and v — 5 m/s) is comparable to the 
j3 decay rate for the very first quantum states, and is more than ten times 
larger for state number bigger than 10. We conclude that specific efforts have 
to be undertaken to set the verticality at an accuracy better than 10~ 5 . This 
will be challenging, because the height of the side mirror is as low as ~ 1 
mm, but such a precision seems possible. 

3.3 Losses due to defects in the edges 

The bottom mirror edges cannot be perfectly flat, we will consider the real- 
istic case where a 50 /zm brink, considered here as a hole, is present for the 
whole boundary. Neutrons can either be lost in this hole, or be reflected by 
the hole in a different quantum state. To estimate the probability of leaving 
the initial state \N) during the collision at the corner, we assume a free fall 
evolution (without bottom mirror) during the classical time: 

2 x 50 fim 

tlree fall = = } = 20 /ZS, (20) 

5 m/s 
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Unstability of quantum levels due to non-verticality of the trap wall 
1 E= ^ 




Figure 4: Escape probability from quantum states while colliding at the side 
wall as a function of the quantum state number, for different values of the 
angle a between side wall and vertical axis. 



the effective size of the hole is twice the geometrical size due to reflection on 
the vertical mirror. Since the propagator of free fall evolution is well known, 
the wave function at time tf ree f a u reads: 

( m \ 1/2 f f -(z+z')gt\ , . . , . . , 

The amplitude for the neutron to be reflected in the same quantum state 
\N) is thus given by the overlap of this evolved wave function with the wave 
function corresponding to \N): 

A N {t) = J ip N (z)i/j(z,t{ ree M \)dz (22) 

The loss probability, 1 — |Aat| 2 , is shown in fig. [5] as a function of the falling 
time, for the first quantum states. This figure shows clearly that the quantum 
states are completely lost after an elapsed time corresponding to 10 fim free 
fall ^2-10 fjLia./ g = 1 ms. It also shows that for smaller times, all quantum 
states have the same probability to be lost, that is, p cornor m 10~ 3 for t free fall . 
As shown in fig. El the related loss rate r^ rner = |.p corner i s ten times larger 
than the (3 decay rate. To reach the ultimate sensitivity, specific efforts 
have to be undertaken to minimize the corner defects. However, even in 
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LOSSES DUE TO FREE FALL IN A HOLE 




t[u-s] 

Figure 5: Probability for losing quantum states due to free fall as a function 
of free fall time, for the first quantum states. 



the pessimistic situation presented above, corner defects do not forbid 10 s 
storage time for trapped quantum levels. 

4 Conclusion 

The GRANIT experiment will measure the transition energies of gravita- 
tionally bound quantum states of neutrons below the kilohertz range, by 
inducing resonant transitions. The commissionning phase of GRANIT will 
start in April, 2008, and the first measurements are expected in 2009. We 
argued that increasing the storage time of trapped neutrons is an essential 
feature of the GRANIT experiment, since the width of the resonance curve 
is inversely proportional to the pulse time. We showed that magnetically 
induced transitions are doable even in flow through mode, and it will be 
possible to resolve the first resonances in the very first stage of the GRANIT 
experiment. The accuracy of the transition energies is limited by the f3 de- 
cay lifetime, and by the shift in energies due to noninertial effects induced by 
Earth rotation. This last subtle effect can in principle be avoided using only 
neutrons travelling in the East- West direction, but this will not be necessary 
since this effect will be dangerous only if we approach the ultimate sensitiv- 
ity. Then the lifetime of trapped quantum levels due to imperfections of the 
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QUANTUM LEVELS LIFETIMES -- GRANIT PRELIMINARY 




v = 5.0 m/s L = 0.3 m 
tan(a) = 1e-05 
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Figure 6: Loss rate of trapped quantum states as a function of the level 
number, due to effects analysed in the text: f3 decay, waviness of the bottom 
mirror, free fall in the corner defect, deviation from verticality of the side 
walls. Also shown the minimum loss rate needed to resolve the state, and 
minimum loss rate needed to be sensitive to Earth rotation bluring of the 
energy levels. 
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trap were estimated, and compared to the ultimate (3 decay storage time. 
We showed that the waviness of the bottom mirror is by no means a prob- 
lem. The main source of losses comes from the side walls of the trap: both 
the deviation from verticality and the corner defects have to be minimized. 
Even in the most pessimistic storage time of 10 s in the trap can be 

reached, corresponding to the relative accuracy 10 -4 of transition energies. 
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